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Abstract 

In the phenomenological description of the nuclear interaction a 
crucial role is traditionally played by the exchange of a scalar 1=0 
meson, the sigma, of mass 500-600 MeV, which however is not seen 
clearly in the particle spectrum and which has a very ambiguous sta- 
tus in QCD. I show that a remarkably simple and reasonably controlled 
combination of ingredients can reproduce the features of this part of 
the nuclear force. The use of chiral perturbation theory calculations 
for two pion exchange supplemented by the Omnes function for pion 
rescattering suffices to reproduce the magnitude and shape of the ex- 
change of a supposed a particle, even though no such particle is present 
in this calculation. I also show how these ingredients can describe the 
contact interaction that enters more modern descriptions of the inter- 
nucleon interaction. 
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When describing QCD to non-physicists, we generally say that it is the 
theory that accounts for nuclear binding. However, in practice our un- 
derstanding of the precise way that QCD leads to nuclear bound states is 
still not good. Nuclear binding is traditionally described by an internu- 
cleon potential which can be parameterized by the exchange of mesons 
The most important exchange producing binding in the central potential is 
a scalar isoscalar meson, the sigma, of mass around 500-600 MeV. While 
other exchanges in the potential are correlated with clear resonances seen 
in the particle spectrum, the sigma is a puzzle. It is not seen in the usual 
way in the spectrum and, after 40 years of debate, does not have a clear 
interpretation in terms of the quarks and gluons of QCD^. It is unfortunate 
that the key ingredient in the signature effect of the strong interactions has 
such an ambiguous status. 

The expectation is that the sigma represents, in some way, the exchange 
of two pions. The quantum numbers certainly are correct for this. So- 
phisticated attempts that construct the potential from scattering data (e.g. 
[3]) have two pions as the lightest intermediate state. However, while phe- 
nomenologically useful, these are not able to answer the question of the 
fundamental nature of the sigma effect. Modern descriptions of the in- 
ternucleon interaction use chiral perturbation theory to calculate two pion 
exchange at low energy [7j. However, these do not produce the sigma effect 
because the chiral amplitudes grow monotonically with the energy. In this 
paper I add a simple and well-motivated addition to the chiral description, 
i.e. the Omnes function describing pion rescattering. We will see that this 
will produce an interaction remarkably close in structure to the exchange 
of a 600 MeV sigma meson. It is clear that there is no resonance in this 
description, yet the needed properties of sigma exchange are reproduced^. 

Recently, there have been successful applications of ideas of effective field 
theory in which the nuclear interaction is treated not by potentials but by 
contact interactions - delta function interactions [71 IH]. At low energy (recall 
that the energy typical of nuclear binding is 10 MeV/nucleon) the result 
of the exchange of a heavy particle can be described by a local interaction. 
Mathematically, this is consistent with the potential description because, 
as the mass m gets large, the Yukawa potential forms a representation of a 
delta function. Physically, this follows from the uncertainty principle, as the 
exchange of a heavy particle has a short range. Nonlocality, to the extent 

discussion of the status of the sigma which is very much in the spirit of the present 
work can be found in A careful recent analysis of tttv scattering describing the sigma 
as a pole on second sheet, quite far from the real axis, is found in 1 

^ Other attempts to describe the nuclear interaction without a sigma are seen in [B| 



1 



it is needed, can then be described by contact derivative interactions. This 
development greatly increases the generality of the description of nuclei, 
as it reduces multiple potentials with different functional forms to a small 
number of constants giving the strengths of the contact interactions. I will 
also show how the use of chiral perturbation theory plus the Omnes function 
can provide a reasonably controlled calculation of the 1 = scalar contact 
interaction. 

In QCD, the low energy interactions of two nucleons obeys an unsub- 
tracted dispersion relation in the different partial wavesjJElEI Specializing 
to the scalar isoscalar channel, we can write the momentum space and co- 
ordinate space interaction as 



Vs{r) = r dfifi e-'^^ psifi) (1) 



Here ps describes the physical intermediate states that occur at energy p in 
the crossed channel. The threshold occurs at twice the pion mass and the 
lowest energy intermediate state is two pions. 

At low energy these spectral functions can be rigorously calculated in 
chiral perturbation theory. The imaginary part of the Feynman diagrams 
describe the physical intermediate states and generate the spectral function 
Ps- For the diagrams of Figs 1 a,b,c these imaginary parts are^l ITU] 

ptiP) = ^ [4ci-^ + csip' - 2ml)] ^^^^^ eip - 2m^) (2) 
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Here ci , C2 , C3 are parameters that describe the NNtttt vertex - these have 
been measured in pion nucleon interactions [Qj [T?H II 1 j . I will address the box 
and crossed box diagrams below. These spectral functions are valid in the 
low energy regime only, and one observes that they grow monotonically with 
the energy. 

However, there is another ingredient which necessarily enters. In the de- 
scription of the TTvr system, unitarity requires the inclusion of vrvr rescattering. 
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Figure 1: Two pion exchange diagrams which arise in chiral perturbation 
theory. 



For a single elastic partial wave, unitarity of the S matrix and analyticity 
require a unique form of the solution, given originally by Omnes[T2]. The 
amplitudes in the elastic region are described by a polynomial in the energy 
times the Omnes function 

U(,) = I ^ (4) 

vr J s s — fj,"^ 

Here 5 is the vrvr scattering phase shift, in our case for the 1=0, J=0 chan- 
nel. Chiral perturbation theory is consistent with this order by order in 
the energy expansion. Following Ref. ^H]) it is known how to match this 
general description to the results of chiral perturbation theory by appropri- 
ately identifying the polynomial. The elastic region in this channel extends 
effectively up to energies of 1000 MeV. 

In practice there has been good success at using the lowest order chiral 
amplitudes, supplemented by the Omnes function. An example close to the 
present problem is 77 — > vrvr in the S wave. Here a lowest order calculation 
supplemented with an Omnes function [141 yields results in close agreement 
with both experiment and with a two loop chiral calculation up to energies 
beyond 700 MeV[IS]. 

I will adopt the Omnes solution matched to the leading order chiral 
result, and will explore possible modifications below. The description of the 
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Figure 2: The left graph gives the pion phaseshifts that are the input into the 
Onmes function, while the right figure shows the real part and the absolute 
square of the Omnes function. 

spectral function then becomes 



The phase shifts can be analyzed in chiral perturbation theory in combi- 
nation with experiment, with the definitive treatment of Colangelo et al 
(CGL)[T6]. Their result for the 1 = 0, J=0 phase shift is shown in Fig 2, 
along with the resulting Omnes functions. Note that there is no sigma res- 
onance in the phase shift near 300 - 600 MeV. A resonance in the elastic 
region is manifest by the phase shift passing through 90 degrees, which cer- 
tainly does not happen near the sigma mass. (If one explores the complex 
plane there is a pole on the second sheet very far from the real axisjl].) In 
producing the Omnes function, I had to extend the phase shifts above the // 
= 850 MeV endpoint of the CGL analysis in order that the principle value 
part of the Omnes function integral be well behaved near the upper end. As 
long as this extension is smooth it has little effect on this calculation. 

With these ingredients, we can display the result for the scalar inter- 
action. In Fig. 3, I show the result for p, along with the individual con- 
tributions of the diagrams of Fig 1. If we had a pure sigma exchange this 
would be delta function at the mass of the a, or a Breit-Wigner shape cor- 
responding to a narrow resonance. One could be forgiven for seeing this 
result as a very broad resonance, even though no resonance exists in the 
formalism. The coordinate space potential is shown in Fig. 3. Also shown 
for comparison is the potential of an infinitely narrow 600 MeV scalar with 
a normalization chosen to match. In practice these are hard to differentiate 



ps{f^) = pfRen{f,) + p's\n{f,)\^ 
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Figure 3: The left figure shows the our results for the spectral function 
p{fi)/lJ, as well as the individual components of diagram 1 a,b,c. The right 
figure shows the coordinate space potential rV{r). There are actually two 
curves in the figure on the right. One is the result of this calculation and 
the second is that of a narrow 600 MeV sigma, with normalization chosen 
to match. The curves cannot be distinguished. 

because the curves are nearly identical. The simple description of Eq. 5 
reproduces closely the spatial variation of the sigma potential. The strength 
of the interaction will be addressed below. 

One can address the robustness of this result by considering possible 
higher order modifications of the basic representation. The NNtttt interac- 
tion has been described by the lowest order chiral Lagrangian. There are 
also energy dependent modifications to these low order results. In partic- 
ular, we expect that there might be form factors depending on the energy. 
To probe this effect, let as modify the tttt interaction by a form factor. 

The choice of n = 1 and m = 800 is shown in Fig. 4. While the relative 
contribution of the two diagrams change (since the latter has the form factor 
squared) and the magnitude is different, the energy variation and spatial 
variation are remarkably similar to the original case. Use of a dipole form 
factor does not change this conclusion. 

We can best address the strength of the interaction by describing the 
magnitude of the contact interaction Gs- For a narrow sigma this would 



(6) 



p{^l)/^l 

-200 
-400 
-600 



0.2 0.4 0.6 0.8 

/X (GeV) 

Figure 4: The spectral function /j(^) / ^ in the presence of higher order energy 
dependent modifications as described in the text 

have the value g^/rn^. There is some uncertainty in the appropriate value 
of Gs, because the preferred fit magnitudes depend somewhat on the cal- 
culational scheme used. However, phenomenological studies of the nuclei 
tend to require Gs = 300 — 450 MeV~^ . The contact interaction is given 
by the strength of the momentum space potential of Eq. 1 evaluated at 

= 0. This is just integral under the integrand shown in Fig. 3. The 
result depend most sensitively on the parameter C3, which is not perfectly 
known. The phenomenological extraction of C3 from ttN data has a large 
error bar, C3 = — 4.7]![} q GeV~^ 121 However, when using an Omnes 
representation, it is likely that this constraint is on the product C3$7(2m,r), in 
which case the value would be C3 = — 3.7^Q g GeV~^. (The other parameter 
choices used were ci = -0.64 GeV-2 and C2 = 3.3 GeV'^, although these 
have only a small impact on the results.) The result for Gs as a function 
of C3 is shown in Fig 4. There is good agreement for the required range of 
magnitudes of Gs for the allowed values of 03 ,? . Here the use of a form 
factor does make a difference. With the form factor described above the 
value of Gs is 40% smaller than without it for a given value of C3. However 
at present understanding this difference may be accounted for by adjusting 
the value of C3. These uncertainties in the appropriate values of C3 and Gs 
keep us from using the magnitude as a precise test of the method. 

In an effective field theory treatment, one keeps pion exchange as an 
explicit degree of freedom while treating the shorter range interactions as 
contact terms. In such a treatment, we should treat the box and crossed 
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Figure 5: The strength of the scalar interaction as a function of the param- 
eter C3. 

box diagrams of Fig 1 dynamicahy, and they should therefore not be in- 
cluded into the contact interaction. For this reason I did not include these 
diagrams in the calculation of the integrand ps{fJ-)/fJ. whose integral gives 
the strength of the contact interaction. However, it should be admitted 
that the different approximations schemes that are used to calculate nuclear 
properties treat the iteration of the one-pion exchange interaction in quite 
different ways. Use of time-dependent perturbation theory, field theoretic 
methods, mean-field methods, Bethe-Salpeter approximations, etc all cap- 
ture different amounts of the box and crossed box diagrams. This is related 
to the serious ambiguity of what is meant by a potential in field theorv[T7]. 
and this likely accounts for the range of fit values of Gs- It seems that 
for the scalar central potential the iteration of the one pion interaction is a 
numerically small compared to the irreducible two pion/sigma contribution, 
for example see ^ and Fig 3.15 of Ref. |2j. Hopefully this ambiguity will 
be cleared up as effective field theory techniques are extended to heavier 
nuclei 18 . 

Chiral perturbation theory plus the Omnes function give a quite simple 
description of the scalar central potential, with a result very similar to the 
exchange of a fictitious sigma particle. However fortunately there is no need 
in such a description to postulate such a scalar particle. This description 
appears to be robust, being qualitatively unchanged by the addition of higher 
order interactions. Besides elucidating a long standing puzzle, these results 
are useful because we have a reasonably solid control over all the ingredients, 
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the chiral amplitudes and the vrvr phase shifts. Questions such as the quark 
mass dependence of nuclear bending may now be addressed in a reasonable 
fashion 19 . The connection of the nuclear interaction to QCD becomes more 
under control. 
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